Spin Zeros and the Origin of Fermi Surface Reconstruction in the Cuprates 
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Two recent quantum oscillation studies find contradictory results concerning the existence of spin 
zeros - zeros of the oscillatory signal induced by Zeeman splitting of the Landau levels. We discuss 
these experiments in light of calculations of the oscillations assuming a spin density wave state. We 
find that the lack of spin zeros reported in one of the experiments is consistent with either hole 
or electron pockets in such a state, if the staggered moment is perpendicular to the external field. 
An analysis for field directions near the planes might be able to differentiate between the two. On 
the other hand, if spin zeros exist as reported in the other experiment, then the staggered moment 
would have to have a substantial longitudinal component. We suggest several experiments to test 
whether this is indeed the case. 

PACS numbers: 74.25. Jb, 72.15. Gd, 75.30. Fv 



Since the pioneering work of Doiron-Leyraud et al.}^ 
quantum oscillation studies have provided new insights 
into the nature of the electronic ground state of the 
cuprates. For overdoped compounds, a large Fermi sur- 
face is observed,'^ consistent with previous photoemission 
studies"^ and paramagnetic band structure calculations. 
For underdoped compounds, this Fermi surface breaks up 
into smaller pockets.H'S'^It is thought that this break up 
is due to a density wave reconstruction of the Fermi sur- 
face, perhaps due to the formation of magnetic stripes.® ^ 
These findings have had a significant impact on our un- 
derstanding of the cuprate phase diagram,!^ though de- 
bates continue about the relation of these small pockets 
to the Fermi 'arcs' observed in zero field photoemission 
studies. 

A key insight into the nature of this underdoped high 
fie ld g round state was recently provided by Sebastian et 
al.^ They looked at the oscillations as a function of 
the angle of the magnetic field relative to the crystallo- 
graphic c-axis. If Zeeman splitting were present, then 
such a rotation study should find certain angles, known 
as spin zeros, where the oscillation amplitude would go 
through zero, with higher angles exhibiting a tt phase 
shift relative to lower ones. Up to an angle near 60°, 
they did not observe this effect. This was confirmed in 
a more extensive study.!-^ As they suggest, such a finding 
would be consist ent with the presence of a spin density 
wave state, '^'21121 a,s we will elaborate on below. Subse- 
quent to this, though, another quantum oscillation study 
was done where zeros from two different frequencies were 
found from fits to the data, with one near 40° and the 
other near 50°.'^ Based on this, the authors suggest a 
paramagnetic ground state instead. 

These results raise a general question about the na- 
ture of the high field ground state, and how much the 
Zeeman splitting is actually reduced if this ground state 
is a spin density wave state. Recently, Ramazashvili,'^ 
following an earlier study by Kabanov and Alexandrov,G5] 
has shown on symmetry grounds^^ that for commensu- 
rate spin density wave ordering, Q ~ {tt,7t), the Zeeman 
splitting should be quenched for both electron and hole 



pockets, as they are centered at the magnetic zone bound- 
ary. This is under the assumption that the magnetic field 
strength is large enough for the spins to be reoriented 
perpendicular to the field. On the other hand, for the in- 
commensurate (stripe) case, only the hole pockets have 
this symmetry. His suggestion was this finding could be 
used to resolve the controversy of whether the observed 
oscillation signal originated from electron or hole pockets. 

In this Rapid Communication, we perform calculations 
for a magnetic stripe in the presence of an external field. 
We verify that in the configuration where the spin mo- 
ments are transverse to the field, the hole pockets do 
not exhibit any Zeeman splitting. On the other hand, 
we find that the Zeeman splitting of the electron pockets 
is strongly reduced compared to that of a paramagnetic 
ground state. Implications of our results in regards to 
the quantum oscillation experiments, and proposals for 
new experiments, will be discussed. 

Our model assumes a linear spin density wave state.'^ 
Without an external field, the secular matrix has di- 
mension N, where the magnetic ordering wavevector is 
Q — {I — 2/N, l)7r. Diagonal elements of this matrix are 
of the form Ck+nQ with n running from 1 to N , where e^ is 
the paramagnetic dispersion. Off-diagonal elements due 
to the stripe order are of the form Vs5n,n±i and VcSn,n±2 
where Vs is the magnetic potential, and Vc its charge 
harmonic (for simplicity, we ignore higher harmonics). 

In the presence of an external field, spin must be ex- 
plicitly taken into account, leading to a secular matrix 
of dimension 2N. Taking the spin quantization axis to 
be along the field direction, the Zeeman splitting enter- 
ing into the diagonal terms would be ^^lbHuz where g 
is the g factor, ^b the Bohr magneton, H the external 
field, and Uz a Pauli matrix. The term Vc would have 
a ctq matrix associated with it instead. Ignoring small 
effects due to spin-orbit coupling except for their renor- 
malization of g, then only the relative orientation of the 
staggered moment and the field is important. Thus, Vs 
becomes Vs (cTx sin (f> + az cos (p) where (j) is the angle be- 
tween the staggered moment and the field. For = 0, 
the staggered moment is longitudinal, and for = 7r/2 



it is transverse. We note that for both the longitudinal 
and transverse cases, the secular matrix decomposes into 
two diagonal blocks with dimension N. The 'down' spin 
block can be obtained from the 'up' spin block by invert- 
ing the sign of H (and Vs in the longitudinal configura- 
tion). For the transverse case, the 'down' spin eigenvalue 
spectrum is equivalent to the 'up' spin one if k is trans- 
lated by Q. For the longitudinal case, with Vc=0, the 
secular matrices are equivalent to the paramagnetic case 
except for a shift of et hy ±^iibH (and thus one has full 
Zeeman splitting). The transverse case is more interest- 
ing; detailed calculations presented below find that the 
Zeeman splitting is strongly suppressed. For each H, we 
determine the chemical potential to obtain the correct 
occupation number (which we take to be 1/iV holes). 

In Fig. la, we show the Fermi surface for N=8 in the 
transverse configuration, with an artificially large value 
of H so as to emphasize the Zeeman splitting. The ef- 
fect of the Zeeman field is to split each of the two hole 
pockets (which have equal area) by shifting their cen- 
ters relative to one another. But the split pockets have 
the same area, as they are related to one another by the 
symmetries of the underlying Hamiltonian. This is asso- 
ciated with the fact that the Zeeman splitting is identi- 
cally zero at the paramagnetic orbit centers, as previously 
derived by Ramazashvili,tSl which he denotes as points 
S = {[2n + 1]/N, l/2)7r. To see this in greater detail, we 
plot in Fig. lb the zero contours of i?„^ — E^^, where 
Ena are the eigenvalues of the secular matrix. For this 
plot, an experimentally relevant value of H is chosen in- 
stead. One can indeed see that the zero contours for each 
band all intersect at the points S. On the other hand, 
these zero contours are not all confined to being near the 
zone boundary ky ^ tt/2. In particular, zero contours 
cut across the electron pockets, and as a consequence, 
even though the two split electron pockets have different 
areas, this difference is much smaller than would occur in 
the longitudinal case. We have not been able to identify 
any symmetry reason for why this occurs. Plots of the 
dispersion along the symmetry axis ky — reveal that for 
the electron band (but not for most of the other bands), 
the Zeeman splitting is very small at the bottom of the 
band, the so-called B pointa^Sat {2n/N,n — l)7r. 

We contrast this with the 7V=10 case shown in Fig. 2. 
The fundamental difference here is that in the paramag- 
netic case, there are now two hole bands, each with the 
same pocket area. Again, the hole pockets are split by 
the Zeeman term, maintaining equal areas. But interest- 
ingly, all bands have a zero contour for i5„^ — Eni equal 
to the ky = tt/2 line. Related to this is the fact that 
the paramagnetic orbit centers of the two hole pockets 
(and related cases with N differing by 4) are not stag- 
gered like for the A^=8 case, but are degenerate (again, 
this degeneracy is connected to the presence of two hole 
bands as compared to the single hole band of the N=8 
case). These points are of the form (4n/iV, l/2)7r. In- 
terestingly, the degeneracy present along this horizontal 
line (ky = n/2) does not occur along any vertical lines. 
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FIG. 1: (Color online) N^8 transverse case with 1/^=0.178 
eV, Vc=0. (a) Fermi surfaces for §77=1000 Tesla (this is 
artificially large so as to emphasize the Zeeman splitting). 
Solid and dotted curves represent the Fermi surfaces from the 
'up' and 'down' spin bands, (b) zero contours of E„f — Eni, 
where E„cr are the eigenvalues of the secular matrix with n 
the band index and a the spin. Here, §77 is taken to be 60 
Tesla. 
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FIG. 2: (Color onhne) iV=10 transverse case with K=0.204 
eV, Vc=0. (a) Fermi surfaces for §77=1000 Tesla (this is 
artificially large so as to emphasize the Zeeman splitting). 
Solid and dotted curves represent the Fermi surfaces from the 
'up' and 'down' spin bands, (b) zero contours of E„f — Eni, 
where E„cr are the eigenvalues of the secular matrix with n 
the band index and a the spin. Here, §77 is taken to be 60 
Tesla. 



We now calculate the Zeeman splitting of the elec- 
tron pocket for the A^=8 and 10 cases, both for the 
longitudinal and transverse configurations. Vs was cho- 
sen to yield an orbit frequency (without Zeeman split- 
ting) for the electron pocket near 530 Tesla as observed 
experimentally.'^ As we use a tight binding fit to band 
theory for ek,we also ran these calculations with a mass 
renormalization, Z^ of 3 to come into better agreement 
with the observed cyclotron mass. This renormalization 
is imposed not only on ek, but on Vg and Vc as well to 
maintain the pocket size. Obviously, the Zeeman effect is 
amplified as Z is increased. Our results are summarized 
in Table I. There, it can be seen that the Zeeman splitting 
is drastically reduced in the transverse case, by factors 
which range from 7.5 to 20 relative to the longitudinal 
case where full Zeeman splitting is observed. We chose 
to do the calculations for a field of 60 Tesla to maximize 



TABLE I: Splitting of the electron orbit frequencies for an 
applied field, |//, of 60 Tesla. A'' is the stripe period, Z the 
mass renormalization, and P the polarization of the spins (L 
for the spins aligned with the field, T for the spins perpen- 
dicular to the field). The quantity 2A.F/ H, where A_F is the 
frequency splitting, is equivalent to gms/m where g is the g- 
factor and rria the 'spin' mass. For N—8, Vs is 0.178 eV; for 
iV=10, Vs is 0.204 eV.I^ 
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10 



P 
L 
T 
L 
T 
L 
T 
L 
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2AF/H 
1.19 
0.06 
3.55 
0.19 
1.30 
0.17 
3.90 
0.52 



the Zeeman effect for physically relevant fields. For the 
A''=8 case in the transverse configuration, we have veri- 
fied that the frequency splitting of the orbits is linear in 
H up to 60 Tesla. 

To compare to experiment, we need to estimate the 
effect the splitting has on the oscillation amplitude. Fol- 
lowing Lifshitz and Kosevich,^^ one finds that the Zee- 
man reduction in the amplitude, Rs, can be written as 
cos{pTrAF/[H cos9]) where p is the harmonic index, AF 
the difference of the Zeeman split orbit areas (expressed 
as a cyclotron frequency) , and 6 the angle of the field rel- 
ative to the c-axis. Comparing to the standard formula 
for i?s, we then equate 2AF/H with gmg/m where nig is 
the 'spin' mass. This is the quantity typically reported 
by experiment. 

In Figs. 3 and 4, we show Rs for the electron pocket 
for the A^=8 and 10 cases, for Z=l and 3, and also for 
both longitudinal and transverse configurations (for the 
fundamental harmonic p—1). For the transverse case, 
spin zeros typically occur at very high angles due to the 
drastically reduced Zeeman splitting in this case. This 
can be contrasted with the longitudinal case, where full 
Zeeman splitting occurs resulting in values comparable 
to a paramagnetic state. We remind that for the hole 
pockets, Rs would be identically 1 for the transverse case. 

We now turn to comparing our results to experiment. 
Sebastian et al^ find no spin zeros up to angles near 
60°, nor any evidence for a tt phase shift of the oscilla- 
tions that would be consistent with passing through such 
a zero. This would definitely be consistent with a hole 
pocket as commented on previously by Ramazashvili,-'^^ 
but could be consistent with electron pockets as well. To 
know for sure would require extending the measurements 
to higher angles, or analyzing higher harmonics .'^Unfor- 
tunately, this is a quantitative issue, since we have shown 
one case in Table I where 2AF/H is as small as 0.06, and 
this could be potentially smaller for different choices of 



0.5- 


L ; \ i| 




T ; : III i 


«' 0^ 


; 1 ;;■; ; 


-0.5- 


""""■■-. ; ; ■ :: 


-1 - 


(a) ■■•-..^ / yp 



1 - 




,, 


0.5- 


L \ 

T \ 


llll 


ai" 0- 




r\m 


-0.5- 

-1 - 


(b) \ 


im 



30 50 

e 



90 



30 



60 



90 



FIG. 3: (Color online) Spin amplitude factor Rs versus angle 
of the field relative to the c-axis with |_ff = 60 Tesla and 
7V=8 for longitudinal (L) and transverse (T) spins, (a) Z—1 
and (b) Z=3, where Z is the mass renormalization. Same 
parameters as Fig. lb. 
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FIG. 4: (Color online) Spin amplitude factor Rs versus angle 
of the field relative to the c-axis with |_ff = 60 Tesla and 
7V=10 for longitudinal (L) and transverse (T) spins, (a) Z—1 
and (b) Z=3, where Z is the mass renormalization. Same 
parameters as Fig. 2b. 



Vs and Ve- 
in contrast, Ramshaw et al!^ were motivated to in- 
clude an Rs factor in their fits in an attempt to describe 
the non-trivial dependence of the overall oscillation am- 
plitude on the field angle. From their fits, they claim ze- 
ros at 50° {gms/m=3.2) for the fundamental frequency 
and 40° (gms/m— 2.1) for a smaller secondary frequency. 
For the fundamental frequency, Rs initially increases as 
6 increases from 0. This would only be consistent with 
a paramagnetic state, or a density wave state of non- 
magnetic origin, such as the d-density wave state.l^Sl It is 
also consistent with the spin density wave scenario if we 
assume that the staggered moment has a large compo- 
nent along the field direction. 

Noting that these two measurements are different 
(Ramshaw et al. measure the c-axis resistance, whereas 
Sebastian et al. measure the change in the skin depth), 
the source of the discrepancy between the two has to do 
with the complex waveform of the oscillations. Even if 
there were only one pocket present, up to four frequen- 
cies could be realized due to (1) bilayer splitting and (2) 
warping of each individual cylinder (leading to extremal 
neck and belly frequencies) . In the analysis of Ref. FTTl 
the authors did not detect any amplitude suppression or 



phase shift that would be consistent with the presence of 
a spin zero. They verified this in more extensive work 
where they also did not detect any such effect in the sec- 
ond harmonic where it would have been more obvious/ 
On the other hand, Ramshaw et a/!^ claim that the pres- 
ence of closely spaced frequencies can mask the presence 
of zeros for each individual frequency. Their analysis in- 
volved a fit including both the primary and a lower sec- 
ondary frequency, though they did not include the higher 
secondary frequency claimed in other measurements.'^^ 
To make further progress would require detailed fits 
including all three frequencies over a wider field range. 
Moreover, these measurements should be supplemented 
by NMR, neutron, and magneto-resistance data to 
look for the orientation of the spins in the high field 
state. So far, little is known about this for underdoped 
YBa2Cu30j, samples in the doping range magnetic oscil- 
lations have been observed. For a lower doping, a field 
induced elastic signal has been observed by neutrons, ^^ 
but the orientation of the spins relative to the field has 
not been determined. Several studies, though, have been 
made on 1/8 doped La2_a;Baa;Cu04, where a magnetic 
stripe state is stable at zero field. What is observed from 
bulk susceptibility studies is a spin flop field of about 6 



Tesla, above which the spins appear to be rotated trans- 
verse to the applied field.'^ Based on this, we feel it quite 
likely for the extremely high fields used in the magnetic 
oscillation studies, the spins are indeed in the transverse 
configuration. But further studies will be necessary to 
see whether this is indeed the case for YBa2Cu3 0j, with 
y ~ 6.5. In the quantum oscillation context, this could 
be tested by comparing angle sweeps in the ac and be 
planes, assuming the spins are locked to a particular axis 
in the plane as occurs for 1/8 doped La2-a;Baa;Cu04 at 
low fields. We note that Sebastian et al. did not find 
any significant chaiiges of their data as a function of the 
in-plane field angle .E! 
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